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The convergence of the Q expansion in the ^^i-^Di mixing angle is investigated through 
the renormahzation condition used to fit the free parameter encountered at NNLO. Although 
the NNLO calculation describes the data well at momenta much less than 50 MeV, a variation 
in the renormahzation condition results in 100% expected uncertainties at momenta as low 
as 100 MeV. This uncertainty can be traced to the presence of large cancelations between 
terms entering at NNLO, which are found to be several times larger than the terms entering 
at NLO. This implies a fundamental lack of predictive power at this order in the Q expansion 
of ei . Given the small empirical size of ei , this lack of predictive power is easily understood 
in terms of simple power counting arguments. 
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I. INTRODUCTION 



Effective field theory (EFT) descriptions of nuclear interactions are intended to provide model indepen- 
dent predictions of observables using systematic expansions based on a separation of scales [^-^. In chiral 
perturbation theory (^PT) (the low energy effective field theory for hadronic interactions arising from QCD) 
this is accomplished by implementing a power counting scheme which treats pionic effects explicitly, while 
heavier mesons are integrated out leaving their effects to be summarized in short distance parameters 
Unlike the cases arising in hadronic physics, the application of to nuclear problems is not straightfor- 
ward, as light scales exist in the problem (such as the scattering lengths) which are formally small from the 
perspective of xPT, but are in fact unnaturally large. 

The original scheme proposed by Weinberg [0 to deal with the unnatural scales implements power count- 
ing at the level of the Lagrangian, and then solves the full Schrodinger equation, effectively iterating the 
interaction to all orders. This scheme, also known as A counting, has enjoyed many phenomenological suc- 
cesses [p|j6 p5| , ^ , |32| . However, due to the nonperturbative nature of this approach, there is no systematic 
power counting at the level of the amplitudes and observables. 

During the past two years, there has been considerable interest in a new way to implement EFT for 
systems of interacting nucleons originally developed by Kaplan, Savage and Wise (KSW) The essence 
of this scheme is to count powers of a low mass scale which is denoted by Q, with 

rriT, ^ Q p^Q (1) 

where a is the NN scattering length, and p is the relative momentum of the nucleon. 

This Q counting scheme is based on an expansion of Q/A, where A is a mass scale characteristic of the 
short distance physics. Since niTr, 1/a and p are formally of the same order the scheme includes pa, rriTra 
and p/rriTr to all orders. The scheme requires that the lowest order NN contact interaction be iterated to 
all orders building in the large scattering lengths. All other interactions including all pionic interactions 
are treated perturbatively. This scheme has been used to calculate numerous observables at next to leading 



order (NLO) with at least modest success [|13|^5 18 19[]. The advantage of the scheme is its rigorous power 
counting at the level of observables which allows for the a priori estimate of the theoretical uncertainties in 
any regime where the Q expansion is converging. This ability to estimate higher order corrections is one of 
the principal motivations for doing calculations with effective field theory, rather than with potential models. 

It is also possible to develop an EFT where the pion is integrated out as heavy. This pion-integrated-out 
theory is essentially equivalent to zero range pseudo-potential theory j|] . It contains an analog of Q counting 
where the light scales are p and 1/a only. To avoid confusion, we will use the phrase, "Q counting," to 
refer to the full theory including explicit pions unless we explicitly indicate otherwise or it is clear from 
the context that we mean the pion integrated out case. Apart from the '^Si-^Di mixing angle, ei, which 
is examined in this paper, calculations done with the pions-integrated-out EFT arc often of comparable 



quality to calculations done to the same order in a theory with explicit pions [g_3 23 2jj31|,^ . Surprisingly, 
for observables related to S'-wave NN scattering, the pion-integrated-out theory fares significantly better 
than the theory with pions explicitly included at momentum less than ^ lOOMeV, even though the p/A 
expansion has a radius of convergence limited by a cut at p — im.^/2 Jl3| , p3| , p4|j3l| ] . 

The fact that the pions-integrated-out theory often does as well, or better than the theory with explicit 
pions raises questions about the convergence of the chiral part of the Q expansion, i.e., the part based on 
the expansion m^/A p3|,p4[. The possible problems with the chiral Q expansion are particularly clear when 
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one examines predictions for observables which are chirally sensitive, in the sense that they are dominated 
by pion exchange physics. Those observables which vanish at some order in the Q counting expansion for 
the pion-integrated-out theory but which are finite at the same order in the theory with cxphcit pions form 
an obvious class of chirally sensitive observables. 

An example of such a chirally sensitive quantity is the "shape function", S, introduced in Ref. It is 
defined as 

S{p^) = pcom - {-I /a + ir,p2) (2) 

so that it removes the first two terms from the effective range expansion and leaves the rest. In Q counting 
at NLO there are two free parameters in S wave scattering which fix a and Ve- The shape function removes 
sensitivity to the fit parameters by subtracting off the first two terms in the effective range expansion. In Q- 
counting with explicit pions, S is order (Q/A)^ due to the pion exchange contributions; formally this can be 
seen in the inclusion ofp/mj^ to all orders. In contrast, S is order (Q/A)"^ in the pion-integrated-out theory. 
Thus, in a regime where Q counting is valid, S is highly sensitive to the pion dynamics. Unfortunately, Q 
counting fails to predict even the sign of S at momenta ~ TOtt- The failure of the theory to describe S at 
NLO suggests that the chiral expansion part of the Q counting may not be working well. 

The conjecture that the chiral part of Q counting is converging poorly has been supported by recent 
calculations by Kaplan and Steele |34[ using a toy potential including explicit one-pion-exchange with a 
realistic scattering length and effective range. In this model, S-wave scattering shows no sign of convergence 
at NNLO. Of course, the model studied was only a toy including only a one-pion-exchange potential. It is 
logically possible that higher-order calculations for the realistic nuclear scattering problem which includes 
non-potential ( "radiation pion" ) effects will behave differently. 

It is important to check the convergence of the chiral part of Q counting. Ideally one should check 
convergence by doing higher-order calculations for the realistic problem including radiation pion effects. 
Unfortunately, at present a full NNLO calculation for s-wave scattering has not been published. However, 
the s-d mixing parameter (which involves the s-wave function at only NLO) has been calculated at NNLO 
by Fleming, Mehan and Stewart (FMS) including all radiation pion effects. The s-d mixing parameter 
should also be a good place to study the chiral part of Q counting as it is chirally sensitive in the sense 
defined above: it is zero at NLO in Q counting for a theory with pions integrated out and nonzero in a 
theory with explicit pions. FMS find extremely good agreement between the NNLO calculation for the usual 
s-d mixing parameter, ei |3^ , and the ei extracted from the data (using the Nijmegen partial wave analysis 
(PWA) [^) for relative momenta up to 150 MeVQ. They conclude that this gives phenomenological support 
for the use of perturbative pions as is done in Q counting. 

In this paper, we investigate the dependence of ei at NNLO on the choice of the renormalization condition 
(i.e., the method used to fit the one undetermined low energy constant). The variations in ei resulting from 
different fits are found to be very large; the variations are of order 100% at p of 100 MeV despite the fact 
that such variations only emerge at order Q/A^ . This sensitivity is a result of cancelations between terms 
entering at NNLO which are much larger than the total value of ei by typically an order of magnitude or 



^In FMS this mixing parameter is denoted by ei, and the mixing parameter in a different parameterization |Q is 
called ei 
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more at momenta above 50 MeV. Moreover one of the large and canceling terms is proportional to a free 
parameter and thus is entirely determined by the fit. This strongly indicates that the large dependence upon 
the renormalization condition is a direct consequence of the lack of predictive power at this order in the Q 
expansion. 



II. S-D MIXING AT NNLO 



In this section the NNLO calculation of ei done by Fleming, Mehan, and Stewart ||3^ is discussed briefly 
and relevant formulae are presented. 

The PDS scheme of KSW is a renormalization scheme in which the Q-counting is manifest. The 
starting point is the effective Lagrangian: 



C = ^Tr (a^S 9^I]t) + ^ Tr(m^S + m,I]t) 
8 4 
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where / is the pion decay constant (approximately 132 MeV in this normalization), uj — rn^/{mu 



(3) 

■m-d) 

is related to the quark condensate and /, ^ = \J e'xjp{i\/2f ■ tt / f) is related to the pion fields in the usual 
fashion in chiral perturbation theory, and the O are nucleon contact interactions. 

The first step is to sum all graphs containing only the contact interaction coupling Co to all orders which 
has the effect of accommodating the unnaturally long scattering length present in NN scattering. As a result, 
the product ap is kept to all orders as required in Q counting. The integral is divergent and is regulated using 
dimensional regularization with PDS. The resulting amplitude is denoted by A^~^\ All other interactions 
are included perturbatively. 

Conventionally p% the S matrix in the spin triplet channel is parameterized in the form 



S = 



cos(2ei) e*(2'^o) 
sin(2ei)e*('^"+'^^ 



sin(2ei)e*(*n+*2) 
cos(2ei)e'(2<5o) 



(4) 



where the above is taken to be the definition of ei. 

At NLO ^0 depends on two parameters which are fit to the low energy scattering data, while ei and 82 
are both predictions. At NNLO {i.e. order (Q/A)^), however, the 6*2° coupling constant, which summarizes 
short distance s-d mixing also enters. Thus at NNLO, there is one new parameter which affects ei apart 
from those fit at NLO. 

FMS derived an expression for e at order (Q/A)^: 
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where A'^ ^Ms the leading order contribution to the S wave scattering amplitude given by: 

47r 1 47r 



7 ■ 



(8) 



M 7 + ' ' MCo 

This amplitude contains a pole a.t p = ij representing the deuteron bound state. Numerically, a binding 
energy of 2.22McV results in 7 = VMB = 45.7McV. 

The s-d mixing amplitude involves one exchange of potential pions, and at order is given by: 

-,2 

(9) 
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(11) 



The amplitude at order Q involves the exchange of radiation pions, as well as potential pions. The resulting 
s-d mixing amplitude at order is given by: 

^2 
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where y{a) and Z{a) are given by: 
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Ci, (^2, and are fXR independent dimcnsionlcss combinations of coupling constants: 
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The coefficients Ci and C2 are fit from the NLO s-wave phase shifts. There is a certain amount of art in 
the fitting of these coefficients. Numerically Mehan and Stewart find (i = 0.2345 and (2 — —0.1038 in ffis 
to the Nijmegen partial wave analysis |39| (PWA) heavily weighted toward low momenta |^^. Note that 
although the preceding amplitudes are order (Q/A)^, from Eq. (|^) they give rise to a contribution to the 
mixing angle, ei at order (Q/A)^. 



III. SENSITIVITY AT NNLO 

The calculation of the '^Si - '^Di mixing angle, ei, contains a single free parameter, C|°, at NNLO (apart 
from those fit for s-wave scattering at NLO). This parameter contains all information about short distance 
mixing entering at this order, and appears in Eqs. (|l^) and (|l5|) as (^q. The linear combinations of Eq. ( p^ ) 
explicitly absorb all /Ir dependence into the undetermined short distance coupling constants. In essence, the 
fits of Ci,2,3 are the renormalization conditions relating the parameters in the theory to physical quantities. 
In principle, one has a certain amount of freedom in choosing precisely which physical quantity is fit. Of 
course, to the extent the Q expansion is well converged to the correct result, one expects that the dependence 
of the final prediction on the particular choice of renormalization condition (i.e., the particular quantity to 
which one fits) will be confined to higher-order terms which are numerically small. Thus, one way to test the 
sensitivity to higher order terms is by testing the sensitivity to imposing different renormalization conditions. 

As noted above, C1.2 can be fit from S wave scattering data at NLO. Thus, the NNLO calculation of ei 
depends on a single free parameter associated with the s-d mixing counterterm. The method used in |3^ 
to fit ^6 is via a global fit to the Nijmegen partial wave analysis of the scattering data, weighted towards low 
momenta; this results in — 0.385. The result of this fit for momenta up to 150 MeV is given in Fig. |^. 
FMS conclude that the high quality of this fit supports the use of perturbative pions as is required in the Q 
counting scheme. 
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3.0 




p( MeV) 

FIG. 1. Results of Q-counting calculations of ei. The solid line is the multi-energy Nijmegan partial wave analysis. 
The short dashed curve is the leading order prediction in Q-counting. The long dashed curve is the full NNLO 
calculation. 

In this paper we investigate the sensitivity of the predicted s-d mixing angle to the choice of renormahzation 
condition in order to assess the contribution of higher-order terms. We exploit the fact that there is only one 
free parameter at this order. Thus any fitting procedure must be equivalent to fixing the predicted value of 
the mixing angle at some particular fixed value of momentum. For example, although the fit of FMS was a 
weighted global fit, the predictions would have been the same had the renormahzation condition for (^e been 
that the NNLO prediction for ei agree with the Nijmegen PWA value at a momentum of approximately 111 
MeV. Thus an efficient way to explore different valid choices for renormahzation conditions is to construct 
fits which ensure that the theory agrees with the data at different fixed choices of momentum p^. 

Assuming even modest convergence of the Q expansion (A > SOOMeV) , it should be reasonable to fit 
at any < 200 MeV as we are working at NNLO. By varying from — 200MeV, a continuous set of 
presumably valid renormahzation conditions is obtained which can be used to estimate the size of theoretical 
uncertainties present at this order due to higher-order terms. In this manner, the range resulting from the 
variation of the pn can be used as a gauge for the predictive power of the theory. The results of this variation 
are displayed in Fig. ^. Figure ^is an indication that the Q expansion is not very predictive at this order. By 
100 MeV, 100% differences are present between the fits at 20 MeV and 200 MeV. By 200 MeV the differences 
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reach 200%. 



30 




p (MeV) 

FIG. 2. Results of the variation of the renormalization condition. The bold line is the multi-energy Nijmegan 
partial wave analysis. The thin lines are the predictions resulting from fits at pa = 20 — 200 MeV, in 20 MeV 
intervals. 

The great sensitivity to the choice of renormalization condition casts doubt on the usefulness of perturba- 
tive pions at least for this quantity at this order. Another indication of a general problem with the approach 
is the overall shape of the curve. It is easy to see that for all p ^ Pr, the EFT prediction overshoots the 
data as it climbs essentially quadratically, failing entirely to predict the slow linear ascent of the data. The 
difficulty in reproducing the detailed shape is not immediately apparent in the work of FMS. By considering 
only comparatively small momenta (150 MeV or less) and choosing a relatively high renormalization condi- 
tion, Pr w 110 — 120 MeV, FMS have left very little room for this serious problem to manifest itself. Note, 
however, it is inconsistent with the power counting rules of Q counting to restrict p^ to momenta below 150 
MeV. The counting rules wore derived for p ~ rriTr. It is extremely implausible that momenta just above 150 
MeV cannot be considered of order m^. 
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IV. CANCELATIONS AT NNLO 



The great sensitivity to the choice of renormahzation condition is in some sense paradoxical. The values 
of ^6 for our fits which cross the data vary only from .359 (for pn of 20 MeV) to .412 (for pji of 200 MeV). 
The apparent paradox is in explaining how such modest variations in the fit parameter can yield such large 
variations in the predicted values. This is particularly hard to understand in light of the fact that the 
contribution due to is subleading and should have modest effects to the total, regardless of its value. 
The resolution is both quite simple and very troubling: there are cancelations between terms entering at 
NNLO which are much larger than the total value ei itself. These cancelations can be seen by breaking the 

(2) . . . . (2) 

NNLO contribution, e\ , into two pieces: the term proportional to (^q, which we denote by e^Q^p)' '^'^^''^ the 

(2) 

subscript CT indicates counterterm , and the remaining piece which we will denote as e . The expression 

(2) 

^l(CT) ■ 

The expression for e^^^ can be extracted from Eqs. (0),(|o|), and (^2|) with Cg set equal to zero. 

In Fig. H we show the ei mixing angle extracted from the Nijmegen PWA, along with the NLO prediction, 
the full NNLO calculation of ei (with fixed using = 100 MeV ), and the two terms which enter at 

(2) -(2) -(2) (21 

NNLO — e^CT) ^1 term and the e^^) term are typically more than an order of magnitude 

larger than the prediction at NLO, the prediction at NNLO, and the result extracted from scattering data. 
Both of these terms are much larger than the leading order terms suggesting serious problems with conver- 
gence. Perhaps more seriously, they are much larger than the overall scale of the final result. Note that 
formally both of these terms are subleading and should each separately give modest corrections to the NLO 
result. Instead, each swamps the NLO, and our prediction depends on delicate cancelations between the 
NNLO terms. The great sensitivity to the choice of renormahzation condition is a result of this requirement 
of very large cancelations. Rather slight changes in the momentum where one fits translate into rather large 
changes in the sum. 
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FIG. 3. ei at NNLO. The NNLO contribution to e\ is divided into two pieces, e^^^' (solid line) and e^c^J (long 
dashed line). The multi-energy Nijmegan partial wave analysis is the dotted line. The NLO and the full NNLO 
calculations are graphed as the short-dashed, and dot-dashed lines, respectively. 

The large scale of tfie cancelations at NNLO suggest that there is no real predictive power at this order. 
The individual quantities which cancel should have corrections of relative order Q/A due to higher-order 
contributions; such corrections can easily be 30%. Since these terms are required to cancel by more than an 
order of magnitude to reach the data, ~ 30% uncertainties in each term translate to ~ 300% uncertainties 
for the total. Thus, there appears to be no predictive power at this order. Indeed, this analysis based on 
the cancelations suggests that lack of predictive power is more severe than one might have concluded on the 
basis of the sensitivity to the renormalization condition found in the previous section. 

One practical consquence of the large cancelations at NNLO and the resultant great sensitivity to the value 
of the parameters is evident in the calculation seen in Fig. ([|). Using the parameters quoted in Ref. |Q we 
produced a plot which differed from the plot in Ref. |33| by a perceptable amount. We believe this small 
difference was a consquence of extremely small rounding errors whose effects were magnified. By adjusting 
the ^6 parameter down from .385 to .3835 we were able to qualitatively reproduce their result. 

Of course, it is logically possible that the large cancelation seen at NNLO is an artifact of the way the 
NNLO contribution was divided into two terms. One can always divide a small number into the difference 
of two unnaturally large numbers; in doing this one can only conclude that there is something wrong with 
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the way the eontribution was divided, and not that there is any problem with underlying theory. However, 
in the present context this appears to be quite unlikely a priori. The division into these two terms was not 
chosen to maximize their individual size, but rather to see the size of the counterterm contribution (which 
by standard power counting arguments contributes at NNLO and should have the characteristic size of a 
typical NNLO contribution). Ultimately what we see is that the scale of the NNLO contribution is much 
larger than the full result and one is depending entirely on the fitting C,q to force the result to be at more or 
less the right size. 

As noted above, it is surprising that the sensitivity of the NNLO prediction of ei to the choice of renor- 
malization condition, although large, is as small as it is given the very large cancelations at NNLO. This 

(2) 

can be understood, in part, in terms of the detailed shapes of the two terms at NNLO. Both the c^ct) ^^'^ 
e'l'-^'' terms are cubic in p for low p, and essentially quadratic for p 100 MeV. Thus the shapes are rather 
similar; if they largely cancel at one momentum they will continue to largely cancel at other momenta. This 
tends to moderate the sensitivity to the choice of renormalization condition. 

The ei extracted from the scattering data is also cubic in p at very low momenta, and with one free 
parameter, it is not surprising that the theory fares moderately well in this regime. However, when fits are 
performed at p^ > 100 MeV, the quality of the fits worsen dramatically - a result of trying to fit an essentially 
quadratic function to a linear one. As the momentum increases , this quadratic function dominates since 
the cancelation between e^^) ^^'^ ^i*'^'' ^"-"^ perfect. Therefore it seems reasonable to conclude that the 
theory is simply not predictive above 100 MeV. 

It is amusing to note that at extremely low momenta <C 50 MeV, all of the fits done with pn as large as 
200 MeV seem to converge to a cubic with approximately the same coefficient at low momentum. Given 
the fact that theory does not appear to be predictive at 200 MeV, we are forced to conclude that this is 
fortutitous. 

V. DISCUSSION 

In the Introduction we raised the question of whether the chiral part of Q counting was converging well. 
In this paper we have seen that ei at NNLO does not appear to be predictive and hence does not support 
the assertion that perturbative pions are useful. It is useful, however, to explore the question of whether 
one would expect the Q counting expansion to be predictive at this order. In fact, even if Q counting were 
converging one might have suspected on the basis of simple power counting arguments that the S-D mixing 
angle would not be predictive at NNLO without actually doing the calculation. 

The leading nonvanishing order contribution to ei is order (Q/A)^. To estimate the expected size of ei one 
needs a characteristic value of Q and A. It is not immediately clear what values one should take for Q or A. 
It is probably sensible to consider ei at a momentum of ~ ttt-tt so that one can take Q ^ m-^. If A 300 MeV 
as suggested by Kaplan, Savage and Wise |l^ the power counting estimate for ei is ^ 25°, which is more 
than an order magnitude larger than the observed value of 2°. Assuming that the Q expansion is sensible, 
this suggests that ei is itself a result of large accidental cancelations. As an unnaturally small observable, 
theoretical difficulties are to be expected, and expansions of ei are expected to converge very slowly if at all. 
For example, corrections to the NNLO are nominally {Q/A)^ corresponding to ^ 6° which is three times the 
observed value, indicating a complete lack of predictive power at this order. 

If Q counting is in fact converging, but ei is just accidently small, one can estimate to what order one 
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must work in order to get an accurate prediction. If one insists on an accuracy of 20%, one estimates 
contributions at their nominal order and A is taken to be 300 MeV, then one has to work to order (Q/A)^, this 
corresponds to next-to-next-to-next-to-next-to-next-to-next-to-next-to-next-to leading order. Thus, even if 
the theory is converging as well as could be expected for A ~ 300 MeV, it is implausible that such calculations 
will ever prove tractable. Of course if A is larger, one need not work to quite as high an order. Nevertheless, 
if A ^ SniTT one needs to work to order (Q/A)^ which also seems to be intractable. 

Naively it seemed that the s-d mixing parameter, ei, should be a good place to test the chiral part of the Q 
expansion, as it is chirally sensitive in the sense defined in the Introduction. However, the unnaturally small 
size of the mixing vitiates its utility as a strong test of perturbative pions. We do note, though, that the very 
fact that the mixing is unnatural in itself does cast some doubt on the use of perturbative pions in general, 
as unnatural results are, by definition, unexpected if the scales in the problem are properly identified. 

In summary, we find that calculations of ci at NNLO do not appear to be predictive as they are highly 
dependent on the fitting procedure used (i.e., on the renormalization condition) and the fitting involves 
delicate cancelations of quantities much larger than the total value. We argue that this lack of predictive 
power is consistent with simple power counting estimates given the small size of ei. Thus, we find that 
calculations of ei to this order do not provide evidence for the use of perturbative pions. 
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